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perturbations. The results also suggest that the topological color code is more robust than the toric 
code, in the parallel magnetic field. 
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I. INTRODUCTION 

Conventional phases of matter such as solids or liquids 
are classically described in the framework of Ginzburg- 
Landau symmetry breaking theory of phase transitions^. 
Over 50 years, it was believed that Landau theory is ca- 
pable of describing any phases of matter by introduc- 
ing an appropriate local order parameter characterizing 
different behaviors of two phases on either side of the 
critical point. Discovery of the fractional quantum Hall 
effect^, introduced an exotic phase of matter carrying a 
new kind of quantum order, topological order, which is 
beyond the scope of the Landau symmetry breaking the- 
ory of phase transitions'^''*. Since then, investigation of 
systems with topological order such as quantum dimer 
models'"'"**, quantum spin models''" and superconduct- 
ing states^^ds became a very active area of research in 
condensed matter physics. Later on, remarkable charac- 
teristics of the topologically ordered system and their ro- 
bust nature motivated the physicists in the field of quan- 
tum information to introduce new models for fault toler- 
ant quantum computation-^°"-^^ . 

The toric code by Kitaev is the first example of error 
correcting quantum models making use of the topologi- 
cally degenerate ground state of the system as a robust 
quantum memory-^''. The ground state degeneracy de- 
pends on the topology of the system and it can not be 
lifted by any local perturbation to the Hamiltonian. Only 
non-trivial topological errors are capable of lifting the 
degenerate ground states of the system^'^. The ground 
states are stabilized by a group of local operators called 
plaquettes and stars making it useful as an error cor- 
recting quantum code. The gapped excitations above 
the topologically-ordered ground states correspond to a 
fascinating class of emerging particles that are neither 
fermions nor bosons. These particles are called anyons. 
Most interestingly, it has been shown that the non-trivial 
braiding of non-Abelian anyons can be used for fault tol- 



erant quantum computation'^*'^^ . 

Topological color codes are a different class of error cor- 
recting quantum codes which allow the implementation 
of the whole Clifford group in a fully topological manner 
in the ground state space. The latter makes the platform 
suitable for quantum computation without resorting to 
braiding of quasi-particle excitations^'^. In contrast to 
the toric code which belongs to the Z2 gauge symmetry 
group, the topological color code is related to the Z2 x Z2 
gauge symmetry group which is a direct consequence of 
the coloring properties of the underlying lattice. 

The robustness of the protected topological order of 
the topological color code is a very important issue when 
it comes to practice. Although the topological degener- 
acy of the ground state can not be lifted by local per- 
turbations, the question arises how much this topologi- 
cal nature is robust against local errors and how large a 
perturbation can be, until the model undergoes a phase 
transition? 

This question has been extensively studied recently 
for the case of the conventional toric code in a uniform 
magnetic field^''"'**. It has been shown that the overal 
phase diagram is very rich containing first- and second- 
order phase transitions, multicriticality^^'^^"'^*, and self- 
duality^''* depending on the field direction. If the transi- 
tion is second order, it is typically in the 3d Ising univer- 
sality class except on a special line in parameter space 
where a more complicated behaviour is detected'^''. Sim- 
ilarly, the phase diagram of a generalized Z3 toric code 
in a field also displays various interesting aspects like a 
3d XY quantum phase transition for the simplest case of 
a single parallel field'*^. 

Such questions have not been studied for the topologi- 
cal color code. It is therefore the main motivation of the 
present work to investigate the phase transition of the 
topological color code in the presence of a magnetic field. 
This is done for a single magnetic field in either x- or z- 
direction on the honeycomb lattice. It is shown that the 
topological color code in a single parallel field is mapped 
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exactly to the Baxter- Wu Hamiltonian'^^ in a transverse 
field on the triangular lattice. To the best of our knowl- 
edge neither of the models is exactly solvable and the 
phase diagrams are unknown. So already the simplest 
case of a single parallel field is fundamentally different 
to the analogous problem for the toric code. In a single 
parallel field, one finds for the latter case a mapping to 
the conventional transverse field Ising model displaying 
a second-order quantum phase transition in the 3d Ising 
universality class^''. 

Here we investigate the topological color code in a 
single parallel field by high-order series expansion tech- 
niques about the low- and high-field limit. Additionally, 
a classical approximation of the Baxter- Wu model in a 
transverse field is performed. Our results clearly reveal 
that the topological phase breaks down to the polarized 
phase through a first-order phase transition. The be- 
haviour of the perturbed topological color code is there- 
fore different to the perturbed toric code. Our results 
suggest that the topological order of the topological color 
code is more robust. 

The outline of the paper is as follows: In Sect. II, the 
topological color codes and their most important prop- 
erties are introduced. In the next section III a magnetic 
field is added to the topological color code. For the case 
of a single parallel field, the mapping to the Baxter- Wu 
model in a transverse field is explicitly given. After- 
wards in Sect. IV, we present the classical solution of the 
Baxter- Wu model in a transverse field already displaying 
a first-order phase transition. The high-order series ex- 
pansions of the low-energy spectrum for both low- and 
high- field limits are the contents of Sect. V. The phase 
diagram of the topological code in a magnetic field is dis- 
cussed in Sect. VI. Finally, Sect. VII is devoted to the 
conclusion. 



II. TOPOLOGICAL COLOR CODE 

Topological color code is a topologically-ordered quan- 
tum stabilizer error correcting code which is defined 
on a certain class of two-dimensional lattices called 2- 
colexes'^^. Local degrees of freedom are spins 1/2 which 
are placed on the vertices of an arbitrary 2D trivalent 
lattice which has 3-colorable faces and can be rapped 
around a compact surface of arbitrary topology of genus 
g. Each vertex is connected to three links with differ- 
ent colors say (red, blue, green) and each link connects 
two faces (plaquettes) of the same color. This colorable 
structure is denoted by TCC {QJ, 2;,?P} where 0, e and 
p correspond to a set of vertices, edges and plaquettes 
respectively. An example of a 2-colex construction on a 
honeycomb lattice is illustrated in Fig. 1. 

The protected subspace C of the model is determined 
by the definition of the stabilizer group il which is gener- 
ated by a set of plaquette operators. For each plaquette 
p, there is a pair of stabilizer operators Xp — (^^gp cr„ 
and Zp = {^ogp cr^ where CTo's are the ordinary Pauli 




FIG. 1: (Color online) A piece of 2-colex on a honeycomb 
lattice placed on a torus. Colored circles at the vertices of 
the lattice represent spins 1/2. A colored plaquette like pi 
is characterized by two closed strings of the other two colors. 
X\ , . . . , Xi are the global strings of the model on a periodic 
structure such as a torus with g = 1. The elementary ex- 
citations are described by open strings like 7 which create 
particles at their end points on the surface of a plaquette (see 
P2) which share the same color with the string. 



operators. A given state G C is left trivially 

invariant under the action of the plaquette operators 
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plaquette operator satisfies the square identity relation, 
(Ap)2 = I = (Zp)2,VAp,Zp e it, leading to eigenval- 
ues ±1. An erroneous state j^E'e) is a state with at least 
one eigenvalue -1 which violates the stabilizing condition 

Ap|*e> = ~|*c> or Zpl^-e) = -|*c)-'^. 

It is possible to construct a Hamiltonian in terms of 
stabilizer operators with topologically degenerate ground 
states where all erroneous states correspond to excited 
states. Considering a trivalent lattice A, the topological 
color code Hamiltonian reads: 



H' 



TCC 



(1) 



where in the following we set J > to be positive and 
we focus on the honeycomb lattice. 



A. String-net operators 

Strings are the generalization of plaquette operators 
which can be either open or closed and defined by the 
product of the Pauli spins of that string. For a generic 
plaquette, such as a blue one (pi) in Fig. 1, two closed 
strings of different colors (red and green) which turn 
around it's boundaries, identify the plaquette. In a simi- 
lar manner, the product of different neighboring plaque- 
tte operators may produce a collection of different bound- 
ary operators. Since all closed strings share an even num- 
ber of vertices, they commute with each other and with 
the plaquette operators. 
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There is also another type of string operators which 
commute with plaquette operators but are not the prod- 
uct of boundary operators. These operators which are 
known as the global loops (fundamental cycles), depend 
on the topology of the manifold on which the system is 
defined and have the character of color. The number 
of these closed loops depends on the genus of the man- 
ifold. For example, a torus with g = 1, has two global 
loops which in contrast to the closed boundaries are non- 
contractible. 

For every homology class of the torus there are two 
closed strings each of one color say green and blue which 
connect the plaquettes of the same color. One can define 
a nontrivial closed strings as'*'"': 



where X\X2X^X\ are defined in Eq. 3 and i,j,k,l = 
0(1) correspond to the global string X appearing (not 
appearing) in the set. These topological global loops can 
take one ground state to another one, and any error of 
this type is not detectable. A generic state of the ground- 
state subspace is then a superposition of the different 
states of the coding space as follows: 



(6) 



Excitations 



cx 



5, 



cz 



(2) 



oei 



where / is the set of spins on the string, /i = 1, 2 stands 
for the homology class of the torus and C characterize 
the color. Following this notation, one can label the non- 
contractible loops of the torus as follows: 



Xi^Si^, X2<^S'[ 



GX 



GX 



BZ 



Z3 ^ 5f ^, 



Za ^ s. 



BZ 



(3) 

For the topological color code, these loops are labeled as 
[Xi . . . X4) and {Zi . . . Z4) and shown in Fig. 1. In the 
next subsection, we will show how these operators man- 
ifest a 16-dimensional subspace for the protected coding 
space C. 



B. Ground state structure 

Since the plaquette operators commute with each other 
and therefore with the full Hamiltonian 1, they share a 
simultaneous eigenstate and the ground state of Eq. 1 is 
readily found by finding the eigenstates of the plaque- 
tte operators. By starting from a polarized state with all 
spins in the up direction lO)**'^' (|2J| is the number of ver- 
tices on the lattice) , the ground state of the system which 
is stabilized by the plaquette operators is constructed. It 
is convenient to denote the un-normalized ground state 
in the following form: 



\^,^)=1[{1 + X,)\0}^ 



(4) 



This state vector is a superposition of all possible closed 
string-nets, a typical feature of the ground states of sys- 
tems with topological order^"-'^''. 

Applying the global loops to the ground state |^gs) 
yields a set of different states which satisfy the stabilizing 
condition. Generically, the set of degenerate protected 
subspace of TCC is defined as: 



C = {\ijkl) : \ijkl) = XlX^Xl Xil^g,)} 



(5) 



Any state j^*) of the system which violates the sta- 
bilizing condition iXp\^>) = |*) and Zpl*) = |vli)), is 
an excited state. For example one can say a state \'^) 
for which Xpl*) = - I*) (Zpl*) = -|*)), has a AT-type 
(Z-type) quasi-particle excitation. Since all of the closed 
strings commute with each other, it is not possible to 
express the excited states of the system in terms of the 
closed strings. The elementary excitations of the TCC 
are described by open strings. Colored strings are open 
if they have endpoints. These endpoints are localized at 
plaquettes which have the same color as the string, see 
also Fig. 1. In particular, a plaquette p2 is an endpoint of 
string 7 if the number of edges of 7 meeting at p2 is odd. 
In terms of string operators, a plaquette p is an endpoint 
of 7 if {Xj,Zp} = or, equivalently, if {Zj,Xp} — 0. 
Thus, open string operators do not commute with those 
plaquette operators at their ends. As a consequence, the 
presence of open strings in the system would violate the 
stabilizing condition producing an elementary excitation: 



1*0 



(7) 



These elementary excitations are Abelian anyons. 
There exist other families of nontrivial topological 
charges which emerge by combining the elementary exci- 
tations. Each single excitation is a boson by itself, as well 
as the combination of two elementary excitations with the 
same color. Furthermore, combination of the excitations 
of different color and type creates particles with semionic 
mutual statistics and two families of fermions'^'. 

Ground states are states where all eigenvalues of pla- 
quette operators are +1 for J > 0. The ground-state en- 
ergy of the system on a lattice with Np plaquettes (corre- 
sponding to 2N sites) is therefore Eq = —NpJ = —2NJ. 
The energy of excited states depends only on the total 
number of —1 eigenvalues of the plaquette operators such 
that each —1 eigenvalue costs the energy 2 J, i.e. the first 
excited state of the system is A^p-fold degenerate (omit- 
ting the topological degeneracy) and it has a total energy 
El = Eq + 2J. Consequently, the energy spectrum is 
equidistant and the elementary excitations correspond- 
ing to —1 eigenvalues are static and non-interacting. 
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III. TOPOLOGICAL COLOR CODE IN 
MAGNETIC FIELD 

Our major objective is to study the robustness of the 
TCC in an external magnetic field. We therefore add 
to the topological color code TJtcc a single parallel field 
in the x-direction. Let us stress that it would be fully 
equivalent to consider the case of a magnetic field in z- 
direction, because only the behaviour of X-type and Z- 
type particles are interchanged. To be specific, we study 
the following Hamiltonian 

H=-jY,{X,^Z^)-h.,Y,al , (8) 

where the first sum runs over all hexagons and the sec- 
ond one over all sites of the honeycomb lattice. In the 
following we set /i^^ > 0. The plaquette operators read 

where the numbers (1, . . . , 6) represent the six sites of the 
hexagon p. 

In the presence of the x field, any trivial closed string 
S'^^ = ^oesC^ pays an excitation energy 2N„hx where 
N„ is the number of vertices of the string. The strings 
stay static and the field perturbation in the Hamiltonian 
8 can be considered as a tension term. In contrast, the 
Zp plaquette operators deform the strings inducing them 
to fluctuate. Hence one can think of them as the kinetic 
term'^'^. As already mentioned in the previous section, 
the excitations are created by open strings. Every open 
strings create two excitations at their ends costing a to- 
tal energy of 4J. In the limit of J » /i^: ~ 0, open 
strings are not energetically favorable and tension is too 
weak to prevent the closed strings from fluctuating. So 
as a result, the ground state of the system is a superpo- 
sition of strongly fluctuating closed string. On the other 
limit, /la; ^ J ~ 0, the strong tension suppresses all of 
the fluctuating strings and the ground state of the sys- 
tem lies over a spin polarized phase in the fleld direction. 
Henceforth, it is physically reasonable to think of a phase 
transition between the two phases in some critical value 
of K/J. 

In order to disclose this phase transition, one has to 
solve Hamiltonian 8 and look at the spectrum of the 
model. Unfortunately, the perturbed TCC is no longer 
exactly solvable. However, the Xp plaquette operators 
commute with the magnetic field and therefore with the 
full Hamiltonian Eq. 8. As a consequence, the eigenvalues 
±1 of all Xp plaquette operators are conserved quanti- 
ties and the system is block-diagonal in these quantum 
numbers. The A"-type excitations are therefore frozen 
with a constant energy gap 2J. It is therefore very rea- 
sonable that the relevant low-energy physics takes place 
in the sub-block where no X-type particles are present, 
i.e. we focus on the sub-block where all eigenvalues of 




FIG. 2: (Color online) Every time a acts on the site o of the 
lattice, it can be considered as the end point of three different 
colored open strings say (red, green, blue) which create three 
Z-type QP on the neighboring plaquettes sharing the same 
color with the strings. The yellow circle denotes the cr^ field 
which is connected to the end points of the red, blue and 
green strings and the three colored circles on the surface of 
the plaquettes are the QP excitations in the topological phase. 
The QPs are created on the vertices of the effective triangular 
lattice. 

Xp plaquette operators are -1-1. The Hamiltonian in this 
sub-block is then given by 

H ^-JNp-J^Zp-K^al . (10) 

P B 

Hence, the magnetic field in x-direction only affects Z- 
type particles. As a consequence, the mutual fractional 
statistics between A"-type and Z-type particles plays no 
role in this sub-block. This fact is the basis for the 
mapping of the perturbed TCC in this sub-block to the 
Baxter- Wu model in a transverse field as discussed in the 
next subsection. 



A. Mapping of the Hamiltonian 

The elementary idea behind the mapping to the 
Baxter- Wu model in a transverse field is to swich from 
the original spin operators to an effective language where 
the eigenvalues ±1 of Zp plaquette operators are the ba- 
sic degrees of freedom. Let us start with the pure TCC 
in the sub-block where all Xp plaquette operators are 
-f 1. As mentionned above, the spectrum of the TCC 
is equidistant and it depends only on the eigenvalues of 
the Zp plaquette operators in a very simple fashion. In 
fact, one can introduce a pseudo-spin 1/2 Pauli matrix 
where the entries on the diagonal exactly correspond 
to the eigenvalues ±1 of the Zp plaquette operators. The 
spectrum of the TCC is then given exactly by 

H^l% = -JNp-jJ2r: , (11) 

i 

where i refers to the center of the plaquettes p. Let us 
note that the centers of the plaquettes build an effective 
triangular lattice (see Fig. 2). We stress that this map- 
ping keeps track of energetic properties. Any informa- 
tion on degeneracies is clearly lost. Indeed, the TCC is 
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mapped to an effective magnetic field which has a unique 
polarized ground state in contrast to the topological or- 
der present for the TCC. 

Next we reformulate the magnetic field term 
—hx X)d '^b terms of pseudo-spin variables. We there- 
fore have to determine how the eigenvalues of Zp pla- 
quette operators are affected by the action of a on 
a certain site D. In fact, the operator just flips the 
three eigenvalues of Zp plaquette operators which are at- 
tached to the site 0. This is a direct consequence of 
the anti-commutation relation between and Zp if is 
part of the hexagon p. In terms of pseudo-spin operators 
the flipping of three eigenvalues corresponds to rf rf 
where the sites 1,2, and 3 build an elementary triangle on 
the effective triangular lattice. Henceforth, the Hamilto- 
nian 10 is mapped to the so-called Baxter- Wu model in 
a transverse field: 

H^-JNp-jY,r:~hxJ2n-Tj-Tk'= , (12) 

i (ijk) 

where the first sum runs over the sites of a triangular 
lattice and the second sum runs over all triangles with 
sites i, j, and k. One should note that the number of sites 
on the effective triangular lattice is equal to the number 
of plaquettes Np on the original honeycomb lattice. The 
magnetic field is therefore mapped to an effective three- 
spin interaction. Let us note that in the dual picture, 
the elementary excitations of the topological color code 
{hx = 0) correspond to simple spin fiips which are static 
and non-interacting hardcore bosons. 



IV. CLASSICAL APPROACH 

Our main goal of this paper is to understand and 
detect the topological phase transition between the 
topologically-ordered phase at low fields and the polar- 
ized phase present at large fields. As a first step to this 
end, we applied a classical approach on the Baxter- Wu 
model in a transverse field. Indeed, this model is much 
better suited for a classical calculation than the TCC 
in a field. The reason is that the topological phase has 
no classical analogue. This is different for Hamiltonian 
12. The TCC is mapped to an effective field term which 
has a polarized ground state being well described by a 
classical configuration where all pseudo-spins point in z- 
direction. Additionally, the ground state of the effec- 
tive three-spin interactions originating from the magnetic 
field term have also a well-defined classical counterpart. 
Indeed, there are exactly four classical ground states for 
J = which are uniquely defined by a three-site tri- 
angular unit cell: | —>■—?—?'), | >), \ < >-<^), and 

I —>■;—;—) where the — (4-) corresponds to the eigen- 
state of the x Pauli operator with -1-1 (—1) eigenvalue. It 
is important to realize that the effect of the effective field 
term on the four classical ground states of the Baxter- Wu 
model is very simple. The field term just rotates every 



spin in essentially the same fashion from the a;-axis in di- 
rection to the z-axis. As a consequence, it is sufficient to 
use a single-site unit cell to the study Baxter- Wu model 
in a transverse field in the classical limit. 

Technically, we replace the vector of pseudo- 
spin Pauli matrices (r^,T^,T^) by classical vectors 
(sin6' cos (/>, sin 6* sin0, cos 6*) with unit length. Conse- 
quently, the classical energy of the Hamiltonian 12 reads: 

^ = -J(l + cos6l) - 2/i^(sin6lcos(/))3 , (13) 

where N — Np is the number of sites on the triangular 
lattice. Note that only two angles (f> and 9 appear in 
this expression because a single-site unit cell has been 
used. In order to get the minimum of the classical energy, 
one finds readily (j) = 0. The remaining expression only 
depends on 9. The minimum is given either by 0tcc = or 
by 9poi depending explicitly on the ratio J/h^. The first 
solution ^tcc corresponds to a fully polarized state in z- 
direction, i.e. it refers to the polarized phase in the dual 
language or the topological phase in the original language 
of the perturbed TCC. This solution is definitely favored 
in the low-field limit <C J. In contrast, the second 
solution 9po\ signals the presence of the symmetry-broken 
phase in the perturbed Baxter- Wu model or the polarized 
phase in the initial problem. It is certainly realized in the 
high-field limit S> J. The phase transition betwen 
these two solution happens at h^/J — 2/^/27 ~ 0.385. 
The classical energy displays a kink and the transition is 
therefore first order. 

The first-order nature of the phase transition in the 
Baxter- Wu model in a transverse field in the classical 
limit is already a hint that the topological phase transi- 
tion of the TCC in a magnetic field is different compared 
to the toric code in a field. The low-energy physics of the 
toric code in a single field corresponds to the transverse 
field Ising model on the square lattice. The latter shows 
a second-order phase transition in the classical limit as 
well as for the full quantum problem. In the next section 
we therefore turn to the original problem of the TCC in 
a magnetic field and we show by high-order series expan- 
sion techniques that the quantum phase transition be- 
tween the topological phase and the polarized is indeed a 
first-order phase transition as suggested by the classical 
approach. 

V. PERTURBATIVE CONTINUOUS UNITARY 
TRANSFORMATION 

The method of continuous unitary transformations 
(CUTs) invented by Wegner *'^""'^ is one of the powerfull 
techniques for diagonalizing or rather block-diagonalizing 
a Hamiltonian. Here we use the perturbative realization 
of the method (pCUTs)^2,43 g^^^y ^^le TCC in an ex- 
ternal magnetic field. The central idea of CUTs is to 
choose a suitable generator for the unitary transforma- 
tion and to (block-) diagonalize the initial Hamiltonian 
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by applying an infinite sequence of successive unitary op- 
erators of particular form. The transformation of the 
Hamiltonian takes therefore place in a continuous fash- 



H(l) = U\l)HU{l) 



(14) 



where I is the flow parameter such that H = H{£ — 0) 
and iJcff = H{£ = oo) is the effective block-diagonal 
Hamiltonian. 

Considering a Hamiltonian of the form H = Hq + xV 
where Hq is diagonal and V is the perturbation, the 
pCUT uses the quasi-particle generator"'^ and it maps 
the Hamiltonian to an effective one which conserves the 
number of quasi-particles. The pCUT is applicable to 
models where the unperturbed part of the Hamiltonian, 
Hq, has an equidistant spectrum. Additionally, the per- 
turbation must be written in the form V = X]!^=-Ar^n 
where T„ increments (or decrements, if n < 0) the num- 
ber of quasi-particles (excitations) by n. Hence, the full 
Hamiltonian reads: 



N 



H = Eo + Q + x 2^ T„ 



(15) 



where Eq is the unperturbed ground state energy, x is 
the perturbation parameter and Q is the hermitian oper- 
ator counting the number of quasi-particles. Most impor- 
tantly, the effective Hamiltonian Hcs conserves the num- 
ber of QP's: [Hf,s,Q] — 0. In the following we demon- 
strate how the low-field and the high-field limit can be 
tackled by pCUTs. 



A. Low-field limit (h^ < J) 

As previously outlined, we aim at a quasi-particle de- 
scription for the low-energy physics of the low-field phase. 
Consequently, one would have to set up a high-order se- 
ries expansion inside the topologically-ordered phase in 
an analogues fashion as it has been done for the toric 
code in a magnetic field'^''^^''^". But this is not neces- 
sary in the present case because one can benefit from the 
mapping of the perturbed TCC to the Baxter- Wu model 
in a transverse field. Indeed, the topological low-field 
phase corresponds to a conventional polarized phase in 
the dual picture. In the following we can therefore ap- 
ply the pCUT method to Eq. 12 in order to study the 
low-field limit of the TCC in a magnetic field. 

The first term in Eq. 12 is an effective field term having 
an equidistant spectrum which is considered as Hq for 
the implementation of the pCUT. Hence single spin flips 
are the elementary excitations above a fully polarized 
ground state. The three-body Ising perturbation changes 
the number of spin flips by n = {±1,±3}. Therefore, 
Hamiltonian 12 can be written as: 



H 
2J 



21 ^^-^ 



■ r_i + T+i 



T. 



+3) 



(16) 



where Q = counts the number of spin flips and 

the Tn operators are defined as 



(ijk) 



(17) 



where bl {b^ ) are creation (annihilation) operators of spin 
fiips on site i. 

The pCUT method maps Eq. 16 to an effective Hamil- 
tonian conserving the number of quasi particles. The 
ground state of the system corresponds to the state with 
zero QP {Q = 0) and the elementary excitation is con- 
tained in the IQP sector {Q = 1). Consequently, we 
can compute the ground-state energy per site Cq, the 
one-particle dispersion uj^^{k) and the one-particle gap 
A'f EE a;'*(fc = 0). For J = 1/2, we have calculated e{f 
and A'^ up to order 10 in the expansion parameter hr^ 



1 1 

500690327 
21432600 



19 4 21436 6 
54 ^ 8505 ^ 
5 74305313819 

' ~ 281302875 
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(18) 



1 - Uhl + 32ht 
18694889252 



134356^ 
81 



893025 



29786981411535707^10 
40507614000 



(19) 



Using standard Pade and dlog Pade approximants as 
well as the bare series of different maximal orders, the 
ground-state energy per site ejf and the one-particle gap 
A'^ are shown in Fig. 3. The bare and extrapolated se- 
ries for the ground-state energy are well converged for 
hx/2J < 0.25. The behaviour of the one-particle gap is 
more complicated. The bare series of the gap has alter- 
nating signs. So only the bare series having a negative 
prefactor for the maximal order (orders 2, 6, and 10) dis- 
play a reasonable behaviour for the one-particle gap (see 
lower panel of Fig. 3). In order to get a clear picture 
about the nature of the phase transition from the be- 
havior of the gap, one need to take a look at the overall 
behavior of the gap in either of the low and high-field lim- 
its. Hence, we postpone any conclusion about the nature 
of the phase transition to the section VI after calculating 
the high-field results. 



B. High-field limit (h^ > J) 

In a large enough magnetic field, the ground state of 
the system is in the polarized phase along the field di- 
rection. To investigate the low-energy excitations of the 
system, the field term of Hamiltonian 8 is the diagonal 
part, Hq, and the topological color code -ffxcc represents 
the perturbation. Again, the overall complexity reduces 
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FIG. 3: (Color online) Ground state energy per site eff (upper 
panel) and the single-particle gap A'^ (lower panel) as a func- 
tion of hx/'iJ in the low-field limit. Solid lines (symbols) rep- 
resent bare series (Pade and dlog Pade approximants) . The 
coupling range between the dashed vertical lines corresponds 
to the location of the first-order phase transition deduced from 
the crossing points of the ground-state energy given in Tab. II. 



strongly due to the fact that the eigenvalues of the Xp 
plaquette operators are conserved quantities. The low- 
energy physics is expected (as already detailed above) in 
the sector where all these eigenvalues are -1-1. In this 
sub-block, the topological color code in a single parallel 
field reduces to Eq. 10. Similarly to the dual formula- 
tion of the low-field limit using the Baxter- Wu model in 
a transverse field, the unperturbed part of the high-field 
limit has a fully polarized ground state and its elementary 
excitations correspond to spin-flip excitations. The per- 
turbation flips all six spins of a hexagon. Regardless how 
many particles already exist on the sites of a hexagon, 
the perturbation Zp can therefore change the number of 
spin flips (QP's) hy n — (0, ±2, ±4, ±6). However, one 
has to keep in mind that the stabilizing condition is 
not violated. In other words, only those processes are 
allowed that flip an 'even' number of spins on a plaque- 
tte, every time the perturbation Zp acts. As a results. 



TABLE I: Action of the T„ operators of the Hamiltonian 20 
on the plaquette states of the lattice. 

T+6 



h\hlblhlblbl\QQQmO) — ^ lllllll) 

T+2 



blblblblb.b. 


000011) — > 


111100) 


blblblbAbe 


000101) — s- 


111010) 


blblbAblh 


001001) — >■ 


110110) 


blbAblblbe 


010001) — >■ 


101110) 


bAblblb% 


100001) — 


011110) 


blblblb,b,bl 


000110) — >■ 


111001) 


blblb,blbA 


001010) — 


110101) 


blb^blb.bl 


010010) — 


101101) 


bAblblb^f^ 


100010) — >■ 


011101) 


blblb,b,blbl 


001100) — 


110011) 


blb^b^.bl 


010100) — 


101011) 


bAblbAbl 


100100) — >■ 


011011) 


bAb-AAbl 


011000) — 


100111) 


bAbAAsbl 


101000) — 


010111) 


b.bAbAsbl 


110000) — > 


001111) 



T-2 = (T+2) 



r-6 = (T+e) 



processes with n = (0, ±4) which involve an odd number 
of spin flips on plaquettes are not allowed. Finally, we 
can write Hamiltonian 10 as follows: 

^ = -N+Q+^^ (r_6 + T_2 + T+2 + T+e) . (20) 

The operator Q = ^l^i counts the number of flipped 
spins. Table I summarizes all the matrix elements of the 
Tn operators in the language of local bosonic operators 
which acts on a plaquette of the lattice. Each state on 
a plaquette is deflned via the conflguration of the six 
spins. The vacuum of the bosonic operators corresponds 
to (1 000000)). Spin flips are represented by I's on the 
sites of the plaquette. 

Bearing the stabilizing condition of Xp in mind, we 
would like to emphasize that although the single spin flip 
channel is the lowest excited state for the diagonal part 
of Hamiltonian 20, it is not the relevant one for the phase 
transition happening between the topological phase and 
the polarized phase. Keeping the Xp operators in their 
ground state enforce that the true lowest excitation of the 
system is the 6-QP sector with an energy gap 6 having all 
six particles located on one hexagon. The next excited 
state is created by acting on the neighboring plaquette 
with an energy gap 8, third one with 10 and so forth. 
Setting hx — 1/2, the ground state energy per site and 
the six-particle gap A^^ = uj^^{k = 0) of the TCC in the 
large-field limit read: 
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,hf 
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144' 
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(22) 



The structure of the T„ operators imply that the six 
QPs have to hop around together on the lattice, from one 
plaquette to the other, and can not disperse individually. 
This situation is similar to the low-field limit where par- 
ticles are created on plaquettes and hop around on an 
effective triangular lattice. In the high-field limit one 
can think of the six spin flip quasi-particles as a single 
compound object which emerges on the plaquettes and 
their kinetics take place again on a triangular lattice. 

Figure 4 illustrates the ground-state energy per site 
Cq^ and the six-particle gap A'^^ as a function of J/2hx- 
The bare and extrapolated ground-state energy curves 
are well converged for J/2hx < 1.4. Similar to the low- 
field gap, the convergence of the bare six-particle gap is 
very poor due to the alternating sign structure of the 
series. Again, we emphasize that by looking at the full 
phase diagram along the whole ranges of field strength, 
one can obtain a better understanding about the nature 
of the phase transition. In the next section we present 
strong evidences for a first-order phase transition. 



VI. PHASE TRANSITION 

To investigate the nature of the phase transition, we 
study the behaviour of the ground-state energy and the 
elementary gap in the whole range of parameters by com- 
bining our results of the low-field and high-field expan- 
sion. Choosing = sinO and J = cos 9, one can join 
low- and high-field quantities in the same plot to get a 
complete picture (Here one has to be careful not to mis- 
take 9 with the one which was defined in section IV. The 
former was used as the angle which minimizes the classi- 
cal energy). 

We start our discussion with the ground-state energy 
per site which is shown in Fig. 5. We find that the 
ground-state energy per site of the low-field and high- 
field limit cross for a certain value 9c giving a kink in 
the ground-state energy per site. Most importantly, the 
location of this kink 9c ~ 0.394 is essentially the same for 
various Fade approximants (see Tab. II) as well as for the 
bare series. The latter can be clearly seen from Figs. 3-4 
due to the fact that the crossing point 9c is located in 
the range of parameters for both limits where even the 
bare series of the ground-state energy is well converged. 
The analysis of the ground-state energy per site points 
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FIG. 4: (Color online) Ground state energy per site eq^ (upper 
panel) and the six-particle gap A*^' (lower panel) as a function 
of J/2hx in the high-field limit. Solid lines (symbols) repre- 
sent bare series (Fade and dlog Fade approximants). The 
coupling range between the dashed vertical lines corresponds 
to the location of the first-order phase transition deduced from 
the crossing points of the ground-state energy given in Tab. II. 



therefore strongly to a first-order phase transition taking 
place at 9c in full accordance with the classical approx- 
imation done for the Baxter- Wu model in a transverse 
field. Surprisingly, the classical phase transition point 
9c ~ 0.367 is very close to 9c- 

This finding is further confirmed by looking at the 
magnetization M = deo/dhx and the ground-state sus- 
ceptibility X = —(P^o/dh^ which can be calculated from 
the ground-state energy via the Feynman-Hellman theo- 
rem. These two quantities are shown in Figs. 6-7. The 
first-order nature manifests itself in a jump in the mag- 
netization at 9c- We stress that the height of the jump 
is essentially the same for the bare series and for various 
Fade approximants confirming quantitatively the loca- 
tion and the nature of the phase transition. 

Finally, we discuss the behaviour of the energy gap A 
as a function of 9 which is displayed in Fig. 8. As already 
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FIG. 5: (Color online) The ground state energy per site eo as 
a function of 6. Solid line corresponds to the bare series of 
order 10. Symbols denote different Pade approximants. 



FIG. 7: (Color online) The ground state susceptibility x of 
the system as a function of 6. Solid line corresponds to the 
bare series of order 10. Symbols denote different Pade ap- 
proximants. 
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FIG. 6: (Color online) The magnetization M of the system 
as a function of 6. Solid line corresponds to the bare series of 
order 10. Symbols denote different Pade approximants. 
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FIG. 8: (Color online) Low- and high-field gap A as a function 
of 6 for different dlogPade approximants. 



discussed, above, the convergence of the bare series is poor 
due to the alternating sign of the series in both limits. 
Additionally, dlogPade extrapolation of both limits show 
no clear tendency in favor of a gap closing signaling a 
second-order phase transition. In contrast, most of the 
extrapolants of both limits cross rather close to 9c as can 
be seen from Tab. II which is in full agreement with a 
first-order phase transition. 

VII. CONCLUSION 

The topological color code is proposed as a platform 
where topological quantum computation can be done in 
the sub-space of the different topological ground states. 



It is therefore interesting to understand how external 
perturbations like a magnetic field affect this topologi- 
cal phase and what kind of topological phase transitions 
do occur. In this work we studied the robustness of the 
topological color code in a single parallel magnetic field 
on the honeycomb lattice. Due to the fact that one sort 
of plaquette operators still commute with the full Hamil- 
tonian, one type of elementary excitations is frozen and 
the relevant low-energy physics takes place in the sub- 
block where only one type of elementary excitations is 
present. 

We have shown that the perturbed color code in this 
sub-block is isospectral to the Baxter- Wu model in a 
transverse field on the triangular lattice. To the best 
of our knowledge this model has not been studied be- 
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TABLE II: Crossing points 0c of the low- and high-field gap 
A (left column) and the ground state energy per site eo (right 
column) signaling a first-order phase transition point. 



A 






eo 


Be 


dlog 


Fade [2, 7] 


0.401 


Fade [2,8] 


0.394 


dlog 


Fade [3, 6] 


0.400 


Fade [3, 7] 


0.395 


dlog 


Fade [4, 5] 


0.375 


Fade [4, 6] 


0.394 


dlog 


Fade [5, 4] 


0.358 


Fade [5,5] 


0.394 


dlog 


Fade [6, 3] 


0.370 


Fade [6, 4] 


0.394 


dlog 


Fade [7, 2] 


0.372 


Fade [7, 3] 


0.395 



fore strongly different to the analogues problem of the 
perturbed toric code. The latter displays a second-order 
quantum phase transition at the critical point h^/J « 
0.328 in the 3d Ising universality class for the case of a 
single parallel field where the low-energy physics of this 
problem is given by the conventional Ising model in a 
transverse field on the square lattice. Our results there- 
fore strongly suggest that the topological color code is 
more robust than the toric code. 
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